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ABSTRACT 


Travel  time  models  are  developed  for  automated  storage/ 
retrieval  (AS/R)  machines.  The  S/R  machine  travels 
simultaneously  horizontally  and  vertically  as  It  moves 
along  a  storage  aisle.  For  randomized  storage  conditions 
expected  travel  times  are  determined  for  both  single 
and  dual  command  cycles.  Alternative  input/output 
(1/0)  locations  are  considered.  Additionally,  various 
dwell  point  strategies  for  the  storage/retrieval 

machine  are  examined.  _ 
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INTRODUCTION 


Over  the  past  years.  Interest  in  automated  storage/ retrieval  (AS/R) 
systems  has  grown  tremendously.  Such  systems  have  Increased  considerably 
In  number  and  probably  will  continue  to  do  so  In  the  future.  The  rapid 
growth  In  the  interest  in  AS/R  systems  can  be  attributed  to  such  benefits 
as  lower  building  and  land  cost,  labor  savings,  reduced  Inventory  levels, 
and  an  Improved  throughput  level,  among  others  [13). 

Much  of  the  early  work  done  to  analyze  AS/R  systems  is  based  on  sim¬ 
ulation.  Bafna  [1]  developed  a  design  package  where  the  optimum  configu¬ 
ration  Is  determined  by  using  simulation  in  conjunction  with  a  search  pro¬ 
cedure.  A  similar  approach  is  presented  by  Koening  [8],  where  the  search 
for  the  "optimum"  configuration  is  limited  to  certain  values  of  the  design 
variables  specified  by  the  user.  Using  simulation  in  designing  such  sys¬ 
tems  has  the  obvious  drawback  of  added  computational  effort,  plus  the 
possibility  of  stopping  with  a  sub-optimum  solution.  The  design  variables 
for  AS/R  systems  typically  can  take  on  a  wide  range  of  values.  Thus, 
simulating  each  feasible  design  long  enough  to  assume  steady-state  behavior 
Is  expensive  computationally  and  undesirable  at  an  early  design  stage  where 
the  user  Is  mainly  Interested  in  obtaining  bench-mark  solutions. 

A  design  package  based  on  analytical  techniques  was  developed  by  White 
and  Bozer  [11).  Based  on  Zollinger's  [12)  cost  model,  the  mathematical 
characteristics  and  relationship  of  the  cost  functions  were  defined  cor¬ 
responding  to  various  elements  in  the  system.  Subsequently,  the  minimum- 
cost  design  was  determined  by  performing  a  Fibonacci  search  over  the  number 
of  aisles  In  the  system. 

Karasawa,  Nakayama,  and  Dohi  [7]  developed  a  cost  model  of  an  AS/RS. 
They  considered  only  single  command  cycles.  Their  total  cost  function  in- 


eluded  terms  for  storage  racks,  S/R  machines,  building,  and  land.  They 
assumed  the  rack  cost  per  opening  to  be  a  linear  function  of  the  height  of 
the  system;  the  S/R  cost  per  machine  was  assumed  to  be  constant;  and  building 
cost  and  land  cost  were  assumed  to  be  linearly  proportional  to  the  floor 
area.  They  solved  the  nonlinear  programming  problem  using  Lagrangian  mul¬ 
tipliers  and  then  obtained  the  best  neighborhood  integer  solution. 

A  number  of  studies  have  been  performed  based  on  simulation;  but  they 
have  emphased  analysis  rather  than  design.  As  examples,  some  studies  have  ad¬ 
dressed  the  throughput/ performance  of  AS/R  systems  under  different  operating 
and/or  storage  policies  assuming  that  the  design  (configuration)  of  the  sys¬ 
tem  is  known.  Among  such  studies  are  those  giver,  by  [2],  [41,  [9],  and  [10]. 

In  this  paper,  closed- form  expressions  are  developed  to  determine  the 
expected  travel  time  associated  with  each  trip  based  on  single  and  dual  com¬ 
mand  cycles.  An  imnediate  application  of  the  travel  time  expressions  is  in 
designing  AS/R  systems  and  in  measuring  the  performance  of  AS/RS  installa¬ 
tions. 


EXISTING  METHODS 

Two  analytical  approaches  have  been  described  for  determining  the  ex¬ 
pected  travel  time.  The  first  approach  is  given  by  Graves,  et  al.  in  [4] 
and  [6].  Their  approach  is  limited  by  the  need  to  assume  the  rack  is  square- 
in-  time.  That  is,  the  dimensions  of  the  rack  and  the  vertical  and  horizontal 
speeds  of  the  storage/ retrieval  (S/R)  machine  are  such  that  the  time  to 
reach  the  row  most  distant  from  the  input/ output  (I/O)  point  equals  the 
time  to  reach  the  most  distant  column,  given  that  the  I/O  point  is  located 
at  the  lower  left-hand  corner  of  the  rack. 
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Empirical  experience  indicates  that  the  optimum  design  for  AS/R  systems 
frequently  is  not  square-ln-time.  For  this  reason  it  Is  desirable  to  deter¬ 
mine  the  expected  travel  time  for  a  rack  that  is  not  necessarily  square-ln- 
time. 

An  alternate  method  of  estimating  the  expected  travel  time  is  given  by 
the  AS/RS  Product  Section  of  The  Material  Handling  Institute,  Inc.  (MHI) 

[13].  The  single  command  expected  travel  time  is  essentially  taken  to  be  equal 
to  twice  the  time  required  to  travel  from  the  I/O  point  to  the  storage  slot 
at  the  center  of  gravity  of  the  rack.  As  can  be  seen  in  Figure  1,  the  single 
command  expected  travel  time  is  equal  to  the  round-trip  travel  time  between 
the  I/O  point  and  point  A.  In  order  to  determine  the  expected  round-trip 
travel  time  for  dual  command  cycles,  it  is  assumed  that  a  storage  is  per¬ 
formed  at  point  A  and  a  retrieval  is  performed  at  point  B,  which  is  assumed 
to  be  located  at  three-fourths  of  the  distance  (horizontally  and  vertically) 
from  the  I/O  point.  As  Illustrated  in  Figure  2,  the  expected  round-trip 
travel  time  for  dual  command  cycles  is  equal  to  the  time  to  travel  from  the 
I/O  point  to  A,  travel  from  A  to  B,  and  then  return  to  the  I/O  point. 

The  MHI  travel  time  model  may  not  provide  an  accurate  representation 
of  the  travel  of  the  S/R  machine.  As  an  example,  when  the  rack  Is  square- 
ln-time  and  randomized  storage  is  used,  the  MHI  model  underestimates  by  25 
percent  the  expected  single  command  travel  time. 

Another  aspect  of  the  MHI  model  that  deserves  examination  concerns  the 
specification  of  the  retrieval  point  Involved  in  the  dual  command  cycle. 

Given  that  the  S/R  machine  has  traveled  to  the  first  opening  for  storage, 
there  is  no  reason  to  assume  that.  In  the  long-run,  the  retrieval  point 
will  be  at  the  assumed  location.  Since  the  retrieval  point  is  assumed  to 


Figure  1:  Determining  the  single  command  expected 
travel  time  using  the  MHI  method. 


Figure  2:  Determining  the  dual  command  expected 
travel  time  using  the  MHI  method. 
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be  three-fourths  of  the  distance  from  the  I/O  point,  the  expected  time  re¬ 
quired  to  return  to  the  I/O  point  will  be  overestimated.  However,  as  noted 
above,  the  expected  time  to  travel  from  the  I/O  point  to  the  storage  location 
will  be  underestimated.  Hence,  it  remains  to  determine  if  the  total  dual 
command  travel  time  is  overestimated  or  underestimated. 

To  facilitate  the  determination  of  the  travel  time  for  an  S/R  machine, 
the  following  notation  is  introduced. 

i  N  *  the  total  number  of  openings  in  the  rack, 
t  j  =  one-way  travel  time  between  the  I/O  point  and  the  i^1 
opening  (tQi  =  tiQ). 

*  one-way  travel  time  between  the  i^  opening  and  the  j**1 
opening  (t^  =  t^.). 

E(SC)  =  expected  single  conmand  round-trip  travel  time. 

E(PC)  =  expected  dual  command  travel  time. 

The  expected  single  command  travel  time  can  be  computed  from  the  following  ex¬ 
pression. 

E<sc>  *  if  j,  2toi  o 

The  expected  dual  command  travel  time  is  given  by 

E(pc>  ■  5p4rx  %  jH1  <‘oi +  +  V  <« 
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Computing  for  each  opening  may  become  time  consuming  for  large  N. 
From  the  double  summation  involved  in  Equation  2,  if  N  =  400,  then  159,600 
terms  must  be  considered. 

A  STATISTICAL  APPROACH 

In  this  section  closed- form  expressions  and  the  methodology  used  to 
develop  the  travel  times  will  be  presented.  The  following  assumptions  are 
made.  • 

1.  The  rack  is  considered  to  be  a  continuous  rectangular  pick  face 
rather  than  an  array  of  discrete  openings. 

2.  The  S/R  machine-  operates  either  on  a  single  or  dual  command 
basis,  i.e.,  multiple  stops  in  the  aisle  are  not  allowed. 

3.  The  rack  length  and  height,  as  well  as  the  S/R  machine  speed 
in  the  horizontal  and  vertical  directions,  are  known. 

4.  The  S/R  machine  travels  simultaneously  in  the  horizontal  and 
vertical  directions. 

5.  In  calculating  the  travel  time,  constant  velocities  are  used 
for  horizontal  and  vertical  travel. 

6.  Randomized  storage  is  used.  That  is,  any  point  within  the 
pick  face  is  equally  likely  to  be  selected  for  storage  or 
retrieval. 

7.  The  I/O  point  is  located  at  the  lower  left-hand  corner  of 
the  pick  face.  Every  trip  originates  and  terminates  at  the 
I/O  point. 

8.  Pick-up  and  deposit  (P/D)  times  associated  with  load  handling 
are  ignored.  The  P/D  times  is  generally  independent  of  the 
rack  shape  and  the  travel  velocity  of  the  S/R  machine. 
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Furthermore,  given  the  load  characteristics,  the  P/D  time 
is  usually  deterministic.  Hence,  it  is  a  straightforward 
matter  to  include  the  P/D  time  after  the  average  travel 
time  has  been  computed. 

As  stated  earlier,  the  pick  face  is  a  continuous  rectangle  with  known 
dimensions  that  can  vary  from  one  application  to  another.  To  "standardize" 
the  pick  face,  let 

s^  »  speed  of  the  S/R  machine  in  the  horizontal  direction, 
sy  =  speed  of  the  S/R  machine  in  the  vertical  direction, 

L  ■  length  of  the  rack,  and 
H  =  height  of  the  rack. 

Now,  let  t^  represent  the  horizontal  travel  time  required  to  go  to  the 
farthest  column  from  the  I/O  station.  Likewise,  let  ty  denote  the 
vertical  travel  time  required  to  go  to  the  farthest  row  (level).  Then, 
by  definition: 


Let, 


and 


T  *  max  {th,  ty} 
b  =  min  {th/T,  ty/T} 


which  Implies  that  0  sb  s  1,  In  subsequent  discussions,  b  is  referred  to 
as  the  "shape  factor".  Note,  if  the  rack  dimensions  and  travel  velocities 
are  such  that  t^  =  ty,  then  b  =  1  and  the  rack  is  said  to  be  "square-in¬ 
time".  Without  loss  of  generality,  assume  that  T  =  th<  That  is,  b  *  t^T. 
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We  will  consider  first  the  single  command  travel  time.  Let  the  storage  (or 
retrieval)  point  be  represented  by  (x,  y)  in  time,  where  0  s.  x  s  1  and 


0  4,  y  s  b.  Travel  from  (0,  0)  to  (x,  y),  say  t  ,  will  be 

Ay 

txy  =  max  (x,  y) 

Now,  let  F(z)  denote  the  probability  that  travel  time  to  (x,  y)  '  less  than 
or  equal  to  z.  That  is, 

F(z)  =  PKt^  s  z) 

Assuming  the  x,  y  coordinates  are  independently  generated 

F(z)  =  Pr(x  s  z)  •  Pr(y  s  z) 


Furthermore,  for  randomized  storage  the  coordinate  locations  are  assumed  to 
be  uniformly  distributed.  Thus, 


and 


Pr(x  s  z)  =  z 


Pr(y  s  z) 


z/b  if  0  <  z  s  b  | 
1  if  b  s  z  £  1  , 


Hence, 


F(z)  = 


z2/b  for  Oizsb 
z  for  b  <  z  s  1 


Therefore, 
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I  2z/b  for  0  s  z  s  b 

f(z)  «  l 

I  1  for  b  <  z  $  1 

Letting  E(SC)  denote  the  expected  travel  time  under  single  command  for  the 
normalized  rack,  then 

,  1  b  1 

j  E (SC)  =  /  z  f(z)dz  =  J  z  f(z)dz  +  /  z  f(z)dz 

z=0  z=0  z=b 

1  b  -  1 

4-  E(SC)  =  Jz^-dz+  jzdz 

z=0  D  z=b 

\  E(SC)  =  i  b2  +  \ 

Therefore, 

E(SC)  =  1  b2  +  1  (3) 

Next,  consider  the  dual  command  cycle.  By  definition,  each  dual  command 
cycle  involves  two  random  locations;  one  representing  the  storage  point,  the 
other  representing  the  retrieval  point.  To  analyze  the  expected  travel  time 
between  the  two  points,  recall  that  any  point  is  represented  as  (x,  y)  in 
time  and  0  s  x  s  1  and  0  s  y  s  b.  Let  tg  be  the  time  required  to  travel 
between  the  storage  and  retrieval  locations.  Also,  let 

F(z)  =  Pr(tg  s  z) 

F(z)  =  Pr{|x1  -  x2|  <  z)  Pr(|y1  -  y2|  *  z) 
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where  (x^,  y^)  and  (xg,  y g )  are  t*1e  tw0  random  points.  First  consider  the 
term  Pr(|x^  -  Xg|  £  2).  Recall,  if  x^nj  are  the  order  statistics 

of  a  sample  Xj,...,  xn,  from  a  population  with  probability  distribution 
function  f(x)  and  cumulative  distribution  function  F(x),  the  difference 
X(n)  -  x(1)  Is  called  the  sample  range  R.  Letting  H(r)  =  Pr(R  s  r),  in  [7] 
it  is  shown  that 


H(r)  =  n  /  f(v)[F(v  +  r)  -  F(v)]n_1  dv 

V=-co 


(4) 


Differentiating  H(r)  yields  the  probability  density  function 


h(r)  *  n(n  -  1)  J  [F(v  +  r)  -  F(v)ln"2  f(v)  f(v  +  r)  dv  (5) 

V=-c° 


Recall  that  0  s  Xj  s  1  and  0  s  Xg  s  1.  With  n  -  2, 


f(x)  = 


for 


F(x)  - 


1 

0 

0 

x 

1 


0  s  x  s  1 
otherwise 

x  <  0 
Osxsl 

X  >  1 


It  should  be  noted  that  f(x)  and  F(x)  are  derived  based  on  the  assumption  that 
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under  randomized  storage  and  x2  are  uniformly  distributed  in  the  interval 
(0,  1).  Since  0  s  R  s  1,  then  0  s  v  £  1  -  r.  Letting  n  =  2,  from  Equation  5 

h(r)  =  z(l  -  r)  (6) 


Letting 


Fx(z)  =  Pr(|x1  -  x2|  *  z) 

then 


FXU)  =  Pr(0  s  R  s  z) 
or 

Fx(z)  =  2  /  (1  -  r)dr 
*  0 

» 

Therefore, 

Fx(z)  =  2z  -  z2 


(7) 


Next  consider  the  term  Pr(Jy^  -  y2j)  *  z,  where  0  s  y .  s  b  for  i  *  1,  2. 
For  this  case 

y  <  0 
Osysb 
y  >  b 

Hence, 
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\  1/b  0  s  y  s  b 

f(y)  *  {  * 

/  0  otherwise 


Substituting  n  =  2  and  rewriting  Equation  5  gives 


b-r 

h(r)  =2  /  f(v)f(v  +  r)dv 
v=0 


b-r  , 

h(r)  =2  /  —k  dv 
v=0  b^ 


>  b(r)  =  K  (b  -  r) 
b^ 


Letting 


Fy(z)  =  Pr(Jy-j  -  y2|)  sz  where  Oszsb 


then 


-7.  for  0  £  z  s  b 
b* 

(8) 

for  b<zil 
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Recall  that  F(z)  =  Pr( tg  s  z)  where  tg  represents  the  "travel  between" 
time.  Based  on  Equations  7  and  8,  F(z)  can  be  written  as: 

F(z)  =  Fx(z)  Fy(z) 

or 

12 

(2z  -  z2)(  l^-%)  for  0  s  z  s  b 
D  b^ 

(2z  -  z2)  for  b  <  z  s  1 

Therefore 

!(2  -  2z)(2z/b  -  z2/b2)  +  (2z  -  z2)(2/b  -  2z/b2)  if  0  s  z  s  b 

. 

2  -  2z  If  b  <  z  si 

Letting  E(T8)  denote  the  expected  travel  time  between  the  two  randomly 
selected  points, 

1  b  1 

E(TB)  =  /  z  f(z)dz  «  /  z  f(z)dz  +  /  z  f(z)dz 
0  O'  b 

or 

E(7S)  *  5  +  l  b2  -  35  b3  (9) 

Let  E(DC)  denote  the  expected  travel  time  for  a  complete  dual  command  cycle 
with  a  normalized  rack.  Then,  by  definition, 

»  13 


E (PC)  =  E(SC)  +  E(TB) 


Hence, 


E(DC)  .§  +  ^b2-^b3 


(10) 


Using  Equation  3  and  10,  the  expected  single  command  and  dual  command  travel 
times  can  be  calculated  for  a  normalized  rack. 

EXAMPLE:  Suppose  the  rack  dimensions  and  the  S/R  machine  speed  is  such  that 
L  -  352  ft. 

H  =  88  ft. 
s^  =  400  fpm 
sy  =  90  fpm 

Using  the  approach  developed  earlier,  we  have 
th  -  L/sh  *=  352/400  *  0.88  mins, 
and 


tv  =  H/sy  =  88/90  =  0.9778  mins. 

Therefore, 

* T  =  max{0.88,  0.9778}  =  0.9778 
and 


b  «  0.88/0.9778  =  0.900 

Hence,  the  "normalized"  rack  is  0.900  mins,  long  in  the  vertical  direction, 
and  1.0  mins,  long  in  the  horizontal  direction.  Using  Equations  3  and  10: 


E(SC)  *  3  b2  +  1  *=1.27  mins. 

E(PC)  =  §  +  \  1)2  ‘  b3  e  l-7140  m1ns- 


14 


To  obtain  the  results  corresponding  to  the  original  rack,  we  "denormal ize" 
the  above  travel  times  to  obtain 

E(SC)  =  .E(SC)  T  =  1.2418  mins. 

E(PC)  *  :E(DC)  T  ,  =  1.6759  mins. 

COMPARISON  WITH  OTHER  METHODS 

As  mentioned  previously,  the  MHI  approach  underestimates  the  expected 
single  command  cycle  time  and  may  overestimate  or  underestimate  the  dual 
command  cycle  time.  First  consider  the  single  command  cycle  time.  For  any 
normalized  rack  with  0  s  b  s  1,  the  MHI  procedure  yields  MHI(SC)  =  0.50  (2)  = 
1.00  minute.  Using  Equation  3  to  determine  E (SC),  the  values  provided  in 
Table  1  are  obtained. 


Table  1.  Comparison  of  Single  Command  Travel  Times  for  a  Normalized  Rack 


a 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0^70 

0.80 

0.90 

H 

E(SC) 

1.003 

1.013 

1.030 

1.053 

1.120 

1.163 

1.213 

H 

1.333 

MHI (SC) 

m 

1.000 

H 

■ 

n 

■ 

*  DIFF 

pm 

0.300 

1.280 

2.910 

5.030 

iomo 

17.560 

1  ■  ■■  , 

21.260 

i 

25000 

It  is  Interesting  to  note  that  the  greatest  difference  occurs  for  b  *  1,  or 
when  the  rack  is  square-in- time.  Also,  the  difference  approaches  zero  as  b 
approaches  zero,  i.e. ,  storage  is  one-dimensional. 

Next  consider  the  dual  command  cycle  time.  For  any  normalized  rack  with 
0  s  b  s  1,  the  MHI  method  yields  MHI(PC)  *  0.75  (21  *  1.50  minutes.  Using 
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Equation  10  to  determine  E(0C),  the  comparison  given  in  Table  2  was  performed. 


Table  2.  Comparison  of  Dual  Command  Travel  Times  for  a  Normalized  Rack 


b 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

m 

E(DC) 

1.333 

1.338 

m 

1.411 

1.454 

1.506 

E 

1.636 

1.714 

1.80 

MHI(flC) 

1.500 

1.500 

1.500 

1.50C 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

1.500 

%  DIFF 

-12.530 

-12.110 

-10.86 

-8.93 

-3.160 

0.400 

4.280 

8.310 

12490 

16.670 

In  Table  2,  a  minus  sign  indicates  over-estimation,  while  a  positive  percent 
difference  indicates  under-estimation  of  the  dual  command  cycle  time.  Notice, 
the  maximum  under-estimation  is  16.67%,  while  the  minimum  is  0.40%  (for  b  * 
0.60). 

From  Equations  3  and  10,  it  can  be  shown  that  a  square-in-time  rack  mini¬ 
mizes  travel  time.  Graves,  et  al.  in  (4)  and  (5]  developed  expressions  for  a 
square-in- time  rack.  Using  a  different  approach,  they  obtained  a  value  of 
7/15  for  the  expected  travel  time  between  two  randomly  selected  points.  Letting 
b  =.  1  in  Equation  9  yields  the  same  value. 

Adjusting  for  Acceleration  and  Deceleration  Rates 

Some  error  is  introduced  when  the  acceleration/deceleration  rate  of  the 
S/R  machine  is  neglected.  Suppose  values  for  s^,  sy  and  the  acceleration/de¬ 
celeration  are  known.  The  total  distance  required  for  the  S/R  machine  to 
start  from  a  stationary  position,  accelerate  to  full  speed  and  immediately 
decelerate  down  to  a  full  stop  can  be  computed.  Suppose  the  distance  is  dj. 


Let  the  time  corresponding  to  dj  be  A  time  units.  The  Impact  of  neglecting 
acceleration/deceleration  is  depicted  in  Figure  3(a).  For  any  distance  less 
than  dp  the  error  is  less  than  A,  and  for  any  distance  greater  than  or 
equal  to  dp  the  error  will  be  A  time  units.  The  total  amount  of  error  in¬ 
troduced  is  represented  by  the  cross-hatched  area.  From  Figure  3(a),  tj  *  th 
(or  ty)  while  tg  is  the  time  required  to  travel  to  the  farthest  column  (or 
row)  from  the  I/O  point.  Including  the  acceleration  and  deceleration  rates. 
One  method  for  reducing  the  amount  of  error  is  to  calculate  b,  E(SC)  and 
E(PC)  as  before  but  to  use  tg  (Instead  of  tj)  in  "denormalizing"  the  rack; 
the  amount  of  error  will  be  reduced  to  the  cross-hatched  area  shown  in  Fig¬ 
ure  3(b).  Assuming  the  acceleration/ deceleration  rates  in  the  horizontal 
and  vertical  directions  are  approximately  equal,  t?  will  not  be  significantly 
different  when  computed  over  L  or  H. 


ALTERNATIVE  CONFIGURATIONS  FOR  THE  I/O  POINT 
In  the  previous  discussion,  it  was  assumed  that  the  I/O  point  is  located 
at  the  lower  left-hand  corner  of  the  rack  and  every  trip  originates  and  ter¬ 
minates  at  the  I/O  point.  In  this  section  the  assumption  is  relaxed;  three 
alternative  configurations  are  analyzed;  and  the  corresponding  expected 
travel  time  expressions  are  developed. 

Input  and  Output  at  Opposite  Ends  of  the  Aisle 

For  the  first  configuration  to  be  considered  assume  that  all  storage 
orders  are  Initiated  at  the  Input  station  while  all  retrieval  orders  are 
terminated  at  the  output  station.  It  is  also  assumed  that  after  each  single 
command  storage,  the  S/R  machine  returns  to  the  Input  point;  and  after  each 
retrieval  (which  may  be  Initiated  at  the  input  point  or  the  output  point) 


travel  j 
time  4 


t 


1 


travel 


the  S/R  machine  travels  to  the  output  station  and  remains  there.  Consequently, 
at  the  time  the  S/R  machine  starts  a  dual  cycle.  It  will  be  at  the  Input  sta¬ 
tion  if  the  previous  trip  was  a  storage;  otherwise  It  will  be  at  the  output 
station.  In  the  first  case,  the  S/R  machine  can  start  the  dual  trip  immediately. 
In  the  latter  case,  before  the  S/R  machine  starts  the  dual  cycle,  it  must  travel 
to  the  input  station.  The  travel  of  the  S/R  machine  is  summarized  in  Table  3. 

Recall,  the  expected  travel  time  for  a  single  command  cycle  when  the  I/O 
is  at  (0,  0)  is 

E(SC)  =  1  b2  +  1 

Dividing  the  above  expression  by  2  gives  the  expected  one-way  travel  time 
E((/),  i.e.,  the  expected  travel  time  from  any  corner  of  the  rack  to  a  randomly 
selected  point  or  vice  versa.  Hence 

E(«0  =  J  b2  +  l 

* 

Also,  recall  the  expected  travel  time  between  two  randomly  selected  points 
is 

E(T8)  ■=  1  ♦  £  b2  -  b3 

Note  that,  E(T8)  will  remain  constant  regardless  of  the  locations  of  the  in¬ 
put  and  output  stations. 

Assume  that  a%  of  all  the  trips  are  single  command  cycles  and  (1  -  a)% 
are  dual  command  cycles.  Furthermore,  since  the  total  number  of  storages 
should  be  equal  to  the  total  number  of  retrievals  in  the  long-run,  assume 
that  50%  of  all  single  command  cycles  are  storages,  while  the  remaining 
50%  are  retrievals.  Also,  assume  that  all  orders  are  statistically 
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Table  3.  S/R  Machine  Travel 


When  the  previous  cycle  was 

To  perform 

S/R  machine  travels 

Single  command  storage 

Single  command  storage 

I+S+I 

Single  command  storage 

Single  command  retrieval 

I+R+0 

Single  command  storage 

Dual  command 

I+S+R+0 

Single  command  retrieval 

Single  command  storage 

0+I+S+-I 

Single  command  retrieval 

Single  command  retrieval 

0-41+0 

Single  command  retrieval 

Dual  command 

O+I+S+R+O 

Dual  command 

Single  command  storage 

0+I+S+I 

Dual  command 

Single  command  retrieval 

O+R+O 

Dual  command 

Dual  command 

O+I+S+R+O 

I:  input  location 
0:  output  location 
S:  storage  location 
R:  retrieval  location 


?0 


independent.  Hence,  the  probability  that  a  given  order  is  a  storage  is  a/2 
while  the  probability  that  it  is  a  retrieval  or  dual  command  order  is  . 

1  -  (<*/2). 


The  above  implies  that,  on  the  average,  a/2%  of  the  trips  will  terminate 
at  the  input  station,  and  1  -  (a/2)%  of  the  trips  will  terminate  at  the  out¬ 
put  station.  If  the  S/R  machine  is  at  the  input  station,  then  the  expected 
travel  time  for  a: 

(1)  storage  is  equal  to  2  E(U) 

(2)  retrieval  is  equal  to  2  E(l/) 

(3)  dual  trip  Is  equal  to  2  E{V)  +  E(TB) 

If  the  S/R  machine  is  at  the  output  station,  then  the  expected  travel 
time  for  a: 

(1)  storage  is  equal  to  2  E(V)  +  K 

(2)  retrieval  is  equal  to  2  E(V) 

(3)  dual  trip  is  equal  to  2  E(t? )  +  E(tb)  +  K 

where  K  is  the  fixed  travel  time  from  the  output  to  the  input  station  (CM). 

Defining  an  operation  as  a  storage  or  a  retrieval,  the  expected  travel 
•  time  per  operation,  say  Ej(T)  will  be: 

El(T)ssf{f  2  E(V)  +f  2  E(l/)  +|  (1  -  a)[2  E(V)  +  E(TB)]| 

+  D  -  f]  {f  U+  2  EU)  +f  2  E(V)  +1  (1  -  a)U  +  2  E(|»  +  E(TB)]} 

which  reduces  to  ! 


21 


E^T)  =f  |  E(l/)(1  +  a)  +1e(T8)(1  -  <*)} 


+  D  _  f D  |  E(C') { 1  +  a)  +  j  E(T8)(1  -  a)  +  ^  K  | 

Hence, 

E-j (T)  =  E(IA)  (1+  ct).  +  1  E(TB)(1  -  a)  +  (1  -  f  )  (11) 

It  Is  instructive  to  note  that  setting  a  =  0  in  the  above  equation  will  yield: 

E^T)  *  E(tO  +^E(7B)  +  \K  (12) 

which  is  intuitively  correct  because  each  dual  trip  involves  two  operations 
(a  storage  and  a  retrieval)  and  total  expected  travel  time  is  E(l/)  +  E(TB) 

+  E(l/)  +  K  =  2  E(l/)  +  E (TB)  +  K. 

Setting  a  =  1  in  Equation  9  gives 

E^T)  *  2  E(K)  (13) 

which  is  intuitively  correct  because  50%  of  all  the  orders  are  storages  and 
the  remaining  50%  are  retrievals;  the  S/R  machine  travels  from  the  output 
station  to  the  input  station  only  if  a  retrieval  is  followed  by  a  storage. 
Note  that  the  S/R  machine  never  travels  from  the  input  station  directly  to 
the  output  station. 

Suppose  the  S/R  machine  is  not  required  to  return  to  the  input  station 
after  a  storage.  Instead,  assume  the  following:  the  S/R  machine  remains  at 
the  point  of  storage,  awaiting  the  next  order.  If  the  next  order  is  a  stor¬ 
age  or  dual  command  order,  the  S/R  machine  returns  to  the  input  station. 
Otherwise,  it  travels  directly  to  the  retrieval  point.  The  travel  of  the  S/R 
machine  is  summarized  in  Table  4. 
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Table  4.  S/R  Machine  Travel 


When  the  previous  cycle  was 

To  perform 

S/R  machine  travels 

Single  command  storage 

Single  command  storage 

X  i 

Single  command  storage 

Single  command  retrieval 

Xj-O^-O 

Single  command  storage 

Dual  comma  '.d 

x,-»-i-*-x2->x3-k) 

Single  command  retrieval 

Single  command  storage 

(KI-+X 

Single  command  retrieval 

Single  command  retrieval 

(KX+O 

Single  command  retrieval 

Dual  command 

0-*- I -*X  i  -*X  2_v0 

Dual  command 

Single  command  storage 

IM-X 

Dual  command 

Single  command  retrieval 

O+X+O 

Dual  command 

Dual  command 

0-*I-*-S+R-»0 

I:  input  location 
0:  output  location 
X:  random  location  in  the  rack 
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Let  I,  0  and  X  denote  the  states  where  the  S/R  machine  is  at  the  input 
station,  output  station  and  at  some  point  within  the  rack,  respectively. 
Then  the  following  transition  matrix  can  be  constructed: 


I 

0 

X 

I 

0 

1  -  - 
2 

a 

1 

0 

0 

>  -I 

a 

2 

X 

0 

1  'I 

a 

1 

For  example,  if  the  S/R  machine  is  in  state  I,  then  it  will  be  in  state  X 
at  the  beginning  of  the  next  trip  only  if  the  current  order  is  a  storage. 

From  Table  4  and  the  above  transition  matrix  it  is  clear  that  the  S/R  machine 
will  never  initiate  a  trip  from  the  input  station;  after  each  storage  the  S/R 
machine  remains  within  the  rack,  otherwise  it  remains  at  the  output  station. 
The  steady-state  probabilities  are  as  follows: 


P 

P 

P 


_  a 
X  '  1 


a 

2 


The  steady-state  probabilities  could  be  obtained  directly  from  Table  4.  How¬ 
ever,  with  complex  operating  procedures,  the  construction  of  a  transition 
matrix  may  simplify  the  analysis. 

Setting  a  =  0  gives  Pq  =  1  and  pj  =  px  =  ^  which  is  intuitively  correct; 
if  only  dual  command  cycles  are  performed,  the  next  trip  will  start  from  the 
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output  station.  On  the  other  hand,  when  a  =  1,  =  Pq  =  ^  and  Pj  =  0  as 

expected,  since  the  number  of  storages  is  assumed  to  be  equal  to  the  number 
of  retrievals  in  the  long-run. 

At  the  start  of  a  cycle,  if  the  S/R  machine  is  at  the  output  station, 
then  the  expected  travel  time  for  a: 

(1)  single  command  storage  is  equal  to  E(i/)  +  H 

(2)  single  command  retrieval  is  equal  to  2  E(t/) 

(3)  dual  command  is  equal  to  2  E(t()  +  E(TB)  +  K. 

If  the  S/R  machine  is  within  the  rack,  then  the  expected  travel  time  for  a: 

(1)  storage  is  equal  to  2  E(tf) 

(2)  retrieval  is  equal  to  E(TB)  +  E(l/) 

(3)  dual  trip  is  equal  to  3  E(^)  +  E(TB) 

Hence,  the  expected  travel  time  per  operation,  E-j(T),  will  be: 

E^T)  =  D  -  f]  Jf  (K  +  E(lO  +  f  2  E(Y)  +  l  (1  -  a) (2  E(l/)  +  E(TB)  +  K)\ 
+  f  Jf  2  E(V)  +  §  (E(TB)  +  E(IO)  +  f  0  -  a)(3  E(lO  +  E(TB))J 


or 

E,(T)-(l-f)  j  j  E(l/)(a  +  2)  +^E(TB)(1  -  a)  +^k| 

♦  f  |  f  E(K)  +  J  (14) 

It  should  be  pointed  out  that  if  all  the  trips  are  performed  on  a  dual  cycle. 
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the  two  strategies  considered  for  the  first  configuration  should  give  iden¬ 
tical  results.  Setting  o  =  0  in  Equation  14  gives 

E^T)  =  E(V)  +  1  E(TB)  +  \  K  (15) 

which  is  identical  to  Equation  12. 

Finally,  it  is  worthwhile  to  compare  the  two  strategies  on  an  expected 
travel  time  basis.  The  second  strategy  is  expected  to  perform  better,  for 
under  the  second  strategy  the  S/R  machine  remains  within  the  rack  and  travels 
to  the  input  station  only  if  the  next  order  is  a  storage.  Assuming  a  *  1, 
from  Equation  14  gives 

E,(T)  =fE(lO  +  +  ]-E(IB)  (16) 

Hence,  from  Equations  IS  and  16,  the  following  ratio,  4,  is  obtained 

Ej(T)  under  strategy  2  ^  E(V)  +  ^  K  +  ^  E(TB)  ^  ^ 

*  ~  Ej(T)  under  strategy  1  2  E[V)  + 1 K 

Assuming  b  =  1,  from  the  equations  for  E(lO  and  E(TB)  It  is  seen  that  E(l /)  = 
o 

2  and  E(TB)  =  7/15.  Furthermore,  for  a  "normalized"  rack,  K  =  1.  Substitu¬ 
ting  these  values  into  Equation  17  gives  $>  =  0.86.  Hence,  for  b  =  1,  using 
the  second  strategy  generates  a  14%  reduction  in  the  expected  travel  time  for 
a  single  command  trip. 
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Input  and  Output  at  the  Same  End  of  the  Aisle,  but  at  Different  Elevations 
The  second  configuration  to  be  considered  is  illustrated  in  Figure  4. 

It  is  assumed  that  the  input  station  is  at  the  lower  left-hand  corner  of  the 
rack  while  the  output  station  is  located  d  time  units  above  the  input  station, 
where  d  <  b.  Furthermore,  it  is  assumed  that  the  vertical  travel  yields  the 
value  of  b. 

Consider  the  single  command  cycle  first.  As  shown  in  Figure  4,  the  rack 
can  be  visualized  as  being  two  separate  racks  (as  indicated  by  the  dashed 
line).  Travel  from  the  input  station  to  a  random  point  in  the  rack  can  be 
expressed  as 

E(V)  =  1  b2  +  \  (18) 

However,  on  going  from  the  random  point  to  the  output  station.  Equation  18 
is  not  appropriate.  The  output  station  may  be  considered  to  be  located, 
at  the  corner  of  racks  A  and  B.  Due  to  symmetry.  Equation  18  holds  as  long 
as  the  station  is  located  at  one  of  the  corners.  Hence,  denoting  the  av¬ 
erage  return  time  from  section  A  as  E^(l/)  and  using  Equation  18  gives 

or 

yr>  -  lLLjp*£. ,  l  (19) 

Similarly, 

Eb(V)  -  £  (20) 

Now,  let  Eq(I/)  denote  expected  travel  time  for  returning  to  the  output  station. 
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Figure  4:  Input  and  Output  at  the  same  end  of  the 
aisle,  but  at  different  elevations. 
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Thus, 


EqOO  *  Pi  <^00)  +  o  -  Pl)<EB(*'))  (21) 


where 


Pj  *  probability  that  the  return  trip  is  Initiated  from  section  A 
Due  to  randomized  storage. 


_  area  of  section  A  _  b  -  d 
P1  "  total  rack  area  ~  b 

Hence,  Equation  21  can  be  given  as 

Simplifying  the  above  expression  gives 

E0(tf)  85  l  b2  +  \  -  |  d(b-d)  (22) 

or 

E000  =  E(V)  -  ^d(b-d) 

Next  consider  the  dual  command  cycle.  By  definition,  travel  time  be¬ 
tween  two  randomly  selected  points  in  the  rack  is  Independent  of  the  location 
of  the  Input  and  output  stations.  Hence,  the  expression  developed  earlier 
for  E(TB)  is  still  valid. 

Assume  that  the  system  operates  according  to  the  first  strategy  defined 
previously.  That  is,  the  S/R  machine  returns  to  the  input  station  after  each 
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storage  and  it  remains  at  the  output  station  after  each  retrieval  or  dual 
command  cycle.  Hence,  if  the  S/R  machine  is  at  the  input  station,  then  the 
expected  travel  time  for  a: 

(1)  storage  is  equal  to  2  E(V) 

(2)  retrieval  is  equal  to  E(V)  +  Eq[V) 

(3)  dual  trip  is  equal  to  E(V)  +  E(TB)  +  Eq(V) 

If  the  S/R  machine  is  at  the  output  station,  then  the  expected  travel 
time  for  a: 

(1)  storage  is  equal  to  d  +  2  E(l/) 

(2)  retrieval  is  equal  to  2  Eg(V) 

(3)  dual  trip  is  equal  to  d  +  E(lO  +  E(TB)  +  Eg(lO 
Thus,  expected  travel  time  per  operation,  E gtT),  will  be: 

E 2(T)  *  f  {  f  2  EOO  +  f  IEOO  +  E0(l/)]  +  \  (1  -  a)[E(U)  +  E(7T3)  +  EQ(lA)]} 

+  (1  ~  f  )  {f  Id  +  2  E(V)J  +|  2.  Eq(|/  )  +  \  (1  -  a) 

.  Id  +  E(tf)  +  E(TB)  +  E0(V)]} 

Simplifying  the  above  expression  yields 

E 2(T)  =  |  {oE(lO  -  \  aE(TB)  +  1  [E(lO  +  E(TB)  +  EQ(i/}]} 

+  (1  -  f  )  {  ?  a  [E(V)  -  E(TB)  +  E0(l/)]  +  \  [E((0  +  E(T8)  +  EgflOl  +  \  d} 

(23) 
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In  Equation  23,  if  o  ■  0  (i.e.,  all  trips  are  dual  cycle  trips),  then, 
as  expected, 

E2(T)  =  \  [E(Pt)  +  E (TB)  +  E0(V)J  +  \  d 

It  is  also  worthwhile  to  compare  the  results  obtained  for  the  second 
configuration  with  those  obtained  for  the  first,  assuming  that  the  same 
strategy  is  used  for  both  configurations.  Recall,  from  Equation  11,  E^(T) 
can  be  expressed  as 

E-j (T)  =  E(U)(1  +  a)  +  \  E(TB)(1  -  a)  +  \  l  (1  -  £  )  (24) 

For  the  purposes  of  the  comparison  let  a  =  0.50,  d  *  0.50b,  and  b  *  1.  By 
definition,  K  =  1.  Thus,  E((A),  E(T$)  and  E0(t{)  can  be  determined  as  follows: 

E(V)  *  l  b2  +  |  «  0.667 

E(TB)  =  \  +  £  bZ  -  Iq  b3  =  0.467 

EQ(I/)  -  l  b2  -  \  bd  +  |  d2  +  1  =  0.542 
Thus,  from  Equation  23: 

E2(T)  *  1.2188  minutes 
and  from  Equation  24:  ' 

E^T)  =  1.4923  minutes 

Hence,  from  a  travel  time  standpoint,  the  second  configuration  performs 
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18.33%  better  than  the  first  configuration.  This  result  was  anticipated,  be¬ 
cause  elevating  the  output  station  will  save  some  travel  time  In  the  vertical 
direction.  Before  selecting  a  configuration  involving  an  elevated  output 
station,  however,  the  costs  associated  with  such  a  design  should  be  considered. 
Input  and  Output  at  the  Same  Elevation,  but  at  a  Mid-Point  in  the  Aisle 

The  third  configuration  alternative  considered  is  based  on  the  I/O  point 
being  located  at  the  center  of  the  rack.  Such  a  configuration  can  be  visual¬ 
ized  as  having  the  delivery  and  take-away  conveyors  running  half-way  into  the 
aisle,  through  a  set  of  rack  openings  located  at  the  mid-level  on  either  side 
of  the  aisle.  It  is  further  assumed  that  vertical  travel  time  e(0,  b),  with 
horizontal  travel  e(0,  1).  Hence,  the  I/O  point  is  assumed  to  be  located  at 
(1/2,  b/2)  for  the  normalized  rack,  where  Os  b  s  1. 

Given  randomized  storage,  the  following  holds: 

Pr(jx  -  s  z) 


2z  if  z  s  1/2 
1  if  z  >  1/2 


Hence, 


Pr(|y  -  ||  s  z) 


f2-  if  z  s  b/2 
1  if  z  >  b/2 


F(z) 


4zZ/b  if  0  s  z  s  b/2 
2z  if  b/2  <  z  s  1/2 

1  if  z  >  1/2 
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Consequently, 


/  8z/b  if  0  s  2  s  b/2 

f(z)  =  <  2  if  b/2  s  z  s  1/2 

\  0  otherwise 

Denoting  expected  travel  time  from  the  center  of  the  rack  to  a  randomly  se¬ 
lected  point  as  EM( V) ,  then: 

1/2  b/2  ?  1/2 

E„(^)  =  /  z  f(z)  dz  «  §  /  z*  dz  +  /  2z  dz 

"  0  0  b/2 

or 

EM(tO  *  ^  b2  +  {  (25) 

Note  that  EM(V)  =  ^E(l/). 

Next  consider  the  strategy  described  for  the  second  configuration.  Since 
the  input  and  output  stations  are  coincident  for  the  third  configuration,  the 
strategy  is  equivalent  to  the  case  where  every  trip  originates  and  terminates 
at  the  I/O  point.  Hence,  the  expected  travel  time  per  operation,  say  E3(T), 
will  be 

E3(T)  =  a{2  Em(I/))  +  (1  -  a) [2  EM(V)  +  E(TB)]  (26) 

It  is  Instructive  to  briefly  compare  the  results  for  the  second  config¬ 
uration  with  the  results  for  the  third  configuration.  Due  to  the  convenient 
location  of  the  I/O  point,  it  is  expected  that  the  third  configuration  will 
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perform  better  from  a  travel  time  standpoint.  However,  the  magnitude  of  the 
Improvement  should  be  determined  in  order  to  economically  justify  higher  con¬ 
veyor  costs  and  the  loss  of  the  rack  openings  required  for  the  conveyor.  Agai 
for  comparison's  sake  assume  o  =  0.50,  d  =  0.50b,  and  b  =  1.  Thus, 

E(y).  =  0.667 
E(T8)  =  0.467 
Eq(V)  =  0.542 
Em((/)  =  0.333 

Hence, 

Eg  CO  =  1*2188  minutes 

and  from  Equation  26 

Ej(T)  =  0.8995  minutes 

which  implies  a  26.2%  reduction  in  expected  travel  time  per  operation. 

Input  and  Output  Elevated  at  the  End  of  the  Aisle 

The  fourth  configuration  alternative  treated  considers  the  situation 
where  the  1/0  station  has  the  location  (0,  d).  As  before,  it  Is  assumed 
that  the  maximum  horizontal  and  vertical  travel  times  are  1.0  and  b,  respec¬ 
tively. 

Recall  the  analysis  of  the  configuration  involving  input  and  output 


stations  at  the  end  of  the  aisle,  but  at  different  elevations.  From  the 
results  obtained,  it  is  straightforward  to  obtain  the  following  expressions 
for  the  expected  travel  times  for  single  command  and  dual  command  cycles: 


(27) 


E (SC)  =  jb2  +  l  -d(b-d) 

E(0C)  =  |  +  1  b2  -  b3  -  d(b  -  d)  (28) 

Comparing  Equations  27  and  29  with  Equations  3  and  10,  elevating  the 
I/O  station  d  time  units  introduces  a  correction  factor  of  d(b  -  d)  in  the 
computation  of  cycle  times. 

DWELL  POINT  STRATEGIES 

The  dwell  point  is  referred  to  as  the  location  of  the  S/R  machine  when 
It  becomes  idle.  In  some  situations,  the  determination  of  the  optimum  dwell 
point  strategy  can  be  treated  as  a  Markov  decision  problem.  In  performing 
such  analyses,  the  travel  time  modeling  approach  used  in  previous  sections 
should  be  applicable. 

The  previous  discussion  focused  on  two  strategies  for  locating  the  S/R 
machine  following  completion  of  storage  and  retrieval  operations: 

A.  Return  to  the  input  station  following  the  completion  of  a 
single  command  storage;  remain  at  the  output  station  fol¬ 
lowing  the  completion  of  either  a  single  command  retrieval 
or  a  dual  command  cycle;  and 

B.  Remain  at  the  storage  location  following  the  completion 
of  a  single  command  storage;  remain  at  the  output  station 
following  the  completion  of  either  a  single  command  re¬ 
trieval  or  a  dual  command  cycle. 

Obviously,  many  other  strategies  could  be  considered.  As  examples  of 
additional  strategies  consider  the  following: 
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C.  Travel  to  a  mid-point  location  in  the  rack  following 
the  completion  of  any  cycle;  and 

D.  Travel  to  the  input  station  following  the  completion 
of  any  cycle. 

Strategy  C  might  be  appropriate  when  the  next  operation  might  be  either  a 
storage  or  a  retrieval;  in  such  a  case,  for  a  given  probability  that  the 
next  operation  will  be  a  storage  an  optimum  dwell  point  location  can  be  de¬ 
termined.  Strategy  D  would  be  appropriate  when  it  is  highly  probable  that 
the  next  operation  will  be  either  a  single  command  storage  or  a  dual  command. 

SUMMARY 

In  summary,  a  number  of  travel  time  models  were  developed  for  automated 
storage/retrieval  machines  under  randomized  storage  conditions.  Expected 
travel  times  were  determined  for  both  single  and  dual  command  cycles.  The 
following  I/O  locations  were  addressed:  I/O  at  (0,  0);  I  at  (0,  0)  and  0 
at  (1,  0);  I  at  (0,  0)  and  0  at  (0,  d);  I/O  at  (0.5,  d);  and  I/O  at  (0,  d). 
Several  dwell  point  strategies  were  considered.  Based  on  the  analyses  per¬ 
formed,  it  is  felt  that  a  number  of  insights  can  be  obtained  concerning 
AS/RS  design  tradeoffs  using  the  travel  time  models. 
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